Observation of Periodic 7r-Phase Shifts in Ferromagnet-Superconductor Multilayers 



A.F. 



in 
o 
o 

(N 
o 

Q 

(N 
(N 

o 
o 

u 

Oh' 



I 

O 

o 



> 
m 

ON 

o 



I 

O 

o 



X 
S3 



V. Shelukhin\ A. Tsukernik^, M. Karpovskii, Y. Bluni\ K.B. Efetov^, 

Volkov'^, T. Champel^, M. Eschrig*, T. Lofwander^, G. Schon"*, and A. Palevski-^ 
^School of Physics and Astronomy, Raymond and Beverly Sackler 
Faculty of Exact Sciences, Tel Amv University, Tel Aviv 69978, Israel 



' University Research Institute for Nanoscience and Nanotechnology, Tel Aviv University, Tel Aviv 69978, Israel 
^ Theoretishe Physik III, Ruhr-Universitdt Bochum, Germany 
Institut fiir Theoretishe Festkorperphysik, Universitdt Karlsuhe, 76128 Karlsruhe, Germany 

(Dated: February 6, 2008) 

We report complementary studies of the critical temperature and the critical current in ferromag- 
net (Ni) - superconductor (Nb) multilayers. The observed oscillatory behavior of both quantities 
upon variation of the thickness of the ferromagnetic layer is found to be in good agreement with 
theory. The length scale of oscillations is identical for both quantities and is set by the magnetic 
length corresponding to an exchange field of 200 meV in Ni. The consistency between the behavior 
of the two quantities provides strong evidence for periodic tt phase shifts in these devices. 
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Devices of superconducting materials as for example 
Josephson junctions 0, have proven exceptionally useful 
in many fields of physics. In bulk conventional supercon- 
ductors, the spin degree of freedom is frozen out by the 
spin-singlet Cooper pair formation. The proximity to a 
ferromagnetic material Q , howeverjOpens up the spin as 
an additional degree of freedom |2|, |j| . Consequently, the 
construction of Josephson junctions of superconductor- 
ferromagnet-superconductor (S-F-S) materials has at- 
tracted considerable attention recently within the emerg- 
ing field of spintronics 

One intriguing effect associated with the new spin de- 
gree of freedom in an S-F-S Josephson junction is the 
thermodynamic stability of a phase difference tt between 
the two superconductors for certain parameter ranges of 
the middle ferromagnetic material 6]. The tt phase is 
a result of a peculiar superconducting proximity effect 
in the ferromagnet (F). The two spin species are split 
in energy by the exchange field E^^i which leads to an 
oscillatory behavior of the proximity induced pair am- 
plitude in the ferromagnet. As a result, properties such 
as the critical temperature Tc and the critical current Ic 
are non-monotonic, oscillating and decaying functions of 
increasing ferromagnetic thickness. For the Josephson 
junction, the free energy loss due to the energy splitting 
can be compensated for by a spontaneous appearance of 
a superconducting phase difference of tt over the junc- 
tion. This additional degree of freedom leads to a series 
of ^ TT and tt ^ transitions, that can be observed as 
zero crossings of Ic and as kinks and minima of Tc with 
varying ferromagnetic layer thickness _7J . 

Previous studies of these effects have been focused 
on ferromagnetic alloys sandwiched between two super- 
conductors, because the corresponding oscillation wave- 
length is quite long and is easily resolved 8, 9J. It is 
important for the development of applications to also un- 
derstand devices including strong ferromagnets as iron, 



cobalt, or nickel, but the short oscillation wavelength 
in these materials is harder to resolve. In previous re- 
ports on devices made of strong ferromagnets either the 
Tc variations [H [ij or the Ic variations IH H [13 
were considered as function of ferromagnetic layer thick- 
ness. However, sample preparation techniques make a 
direct comparison of the different experiments difficult. 
For a successful theoretical understanding a high control 
of material properties is required. 

In this Letter we report studies of both the critical cur- 
rent and the critical temperature variations as function of 
the F layer thickness in S-F-S junctions made of one set 
of materials, namely Nb-Ni-Nb junctions prepared un- 
der identical conditions. We find that both quantities, 
Ic and Tc, vary on the same scale, the magnetic length 
Lm — \/ fiD F I Ecx set by the properties of the ferro- 
magnet only {Dp is the diffusion constant in the ferro- 
magnet). By a detailed comparison of our measurements 
with theory, we find consistent fits for an exchange field 
of Ecx = 200 meV in Ni. 

The samples for the Ic measurements had a 10 x 10 /im^ 
cross sectional area and were fabricated with a standard 
photolithography technique. The process contained three 
stages of lithography: liftoff of the bottom Nb/Cu layers, 
liftoff of the variable thickness Ni layer, and liftoff of the 
top Cu/Nb layers. We fabricated two sets of nine Nb- 
Cu(Au)-Ni-Cu(Au)-Nb junctions with variable Ni thick- 
ness in the range of 35 A to 75 A in steps of 5 A (set I) 
and 4 A (set II). The thickness of each Nb layer is 2000 A, 
while the total thickness of the Cu is 2400±250 A (Au - 
500±50 A). We show the layout of the junctions for the 
Ic measurements in Fig.l 

The Nb films were sputtered using a magnetron gun 
and were covered in situ with a Cu (set I) or a Au (set 
II) layer by thermal evaporation to prevent Nb oxidation. 
The ferromagnet layers of Ni were e-gun evaporated in a 
separate vacuum chamber at a pressure of 2-10~^ torr and 
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FIG. 1: Schematic layout: a) of the junctions for Ic experi- 
ments, b) of planar samples for Tc experiments. 



subsequently covered in situ by Cu or Au. The variation 
of the Ni thickness was achieved by a specially designed 
shutter, which exposed the samples in sequence, so that 
every sample was exposed to the evaporating Ni for ad- 
ditional fragments of time. Because all samples within 
one set were prepared simultaneously, all layer interfaces 
are nominally identical and the only difference between 
the samples is their Ni thickness. The critical current 
was measured by passing a DC current with a small AC 
modulation through the sample. The AC voltage, which 
appeared above the critical DC current, was picked up 
by a lock-in amplifier operated in a transformer mode. 
The measurements were performed in a ^He cryostat in 
the range from 4.2 K down to 1.5 K. 

The samples for the Tc studies were prepared by an in- 
situ evaporation of Ni and Nb layers without photolithog- 
raphy, see Fig. lb. Two sets of structures, each contain- 
ing 16 samples were fabricated. The first set of samples 
contains only a single layer of Nb and was obtained by 
sequential deposition of Ni(30A)-Nb(430A)-Ni(x), where 
Ni(x) was varied from to 37 A in steps of ~2.5 A. 
The second set contains two Nb layers, namely Ni(30A)- 
Nb(430A)-Ni(x)-Nb(430A)-Ni(30A), and was prepared in 
a similar manner. The thicknesses of the bottom and top 
Ni layers, as well as the thicknesses of the Nb layers, were 
chosen such that the bulk Tc of Nb was suppressed. This 
increases the sensitivity of Tc to variations of the thick- 
ness of the center Ni layer. 

For strong ferromagnets like Ni, the magnetic length 
Lm is below 20 A. The first — > tt transition is there- 
fore expected to occur for very thin films. However, Ni 
films thinner than a few tens of Angstroms prepared by 
standard e-gun evaporation is not expected to be homo- 
geneous or to perfectly cover a metallic surface. This 
is a problem for measurements of Ic in S-F-S junctions 
with a very thin F-layer, since uncovered regions short- 
circuit the junction. For thermodynamic measurements. 



such as measurements of Tc, an inhomogeneous coverage 
is less of a problem. At the same time, the effect we ad- 
dress decays exponentially with the layer thickness, and 
it is undesirable to have very thick films, in particular for 
thermodynamics measurements. With this in mind, we 
have chosen Ni-film thicknesses ranging from to 35 A 
for the Tc measurements and thicknesses ranging from 35 
to 75 A for the Ic measurements. 

The results of the critical current measurements are 
shown in Fig. 2. We have also included the set published 
earlier in Ref. Since the sets were prepared sepa- 

rately and are somewhat different in their normal metal 
constituents, they are expected to have different inter- 
face properties reflected as a different transmission co- 
efficient T of the superconductor-ferromagnet interface, 
which enters the calculations of Ic as a prefactor |0| . 
Therefore, we have normalized the values of Ic by 
varying between the sets, as indicated in the legend of 
Fig. 2. We use the Landauer formula to estimate a lower 
limit for T, under the assumption that the entire resis- 
tance of the junction arises from the two S-F interfaces. 
By using a typical resistance value R = 200 fi^l for our 
junctions we get Tmin ~ 0.012. The values of T used for 
all of our sets are indeed larger than the estimated lower 
limit. We emphasize that irrespective of which theory 
our experimental data is compared with, every « 15 A of 
Ni the phase over the Nb-Ni-Nb junction changes from 
to TT (or TT to 0) . This implies that we should expect 
the first pronounced minimum in the variation of Tc at a 
Ni layer thickness of dmin ~ 15 A. 

The theoretical curve for the critical current at 
T=4.2 K normalized by (solid line in Fig. 2) was 
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FIG. 2: The critical Josephson current Ic as function of the 
Ni- layer thickness diMi. The solid line shows the theoretical 
curve with E^x ~ 200 meV used as a fitting parameter. The 
dots are the experimental data. 
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calculated with the theory of Ref. [I3 using the foUow- 
ing parameters: the Fermi velocity vp = 2.8 • 10^ cm/s 
|l7| . a critical temperature of Nb — 8.5 K, the ex- 
change energy E^x = 200 meV, and the mean free path 
^ = 28 A. The magnetic length in terms of these param- 
eters can be determined to be Lj^i = 10 A. Our result for 
the exchange energy is somewhat higher than the results 
obtained from spin-resolved photoemission spectroscopy 
where the splitting between spin-up and spin-down bands 
at the Fermi energy ranged between 2Eex = 200 meV and 
350 meV ^^l- The theoretical predictions for 2Eex are 
typically higher and range from 600 - 850 meV 0| . Our 
value of 2Eex = 400 meV falls in between the above ex- 
perimental and theoretical values, and is consistent with 
the values in Ref. ^3 of Eg^ = 220 meV (corresponding 
to a magnetic length of 8.8 A), obtained from measure- 
ments of Tc in Ni-Nb bi-layers. 

The theoretical prediction of the minima in Fig. 2 were 
obtained with Eq. (19) in Ref. The oscillation period 
for large thicknesses is thus not equal to the thickness 
where the first minimum in Ic is predicted by the theory. 
This is due to the fact that the first minimum occurs at 
a thickness smaller than the mean free path i. Equal 
spacing of the minima only takes place in the regime 
d > i. The theoretical prediction for the first minimum 
using the above fit is dmin « 17 A. 

It was noted recently that the critical current of junc- 
tions with a given thickness of the ferromagnet Py, scat- 
ter considerably We confirm this effect in our Ni 
junctions. We believe that this phenomenon is related to 
the domain structure of the ferromagnet. It was recently 
shown that variations in the domain configuration lead 
to considerable variations in /c, provided that the mag- 
netic flux through typical domains is of the order of the 
flux quantum $0 |20j. By using a magnetization for Ni 
of Ms = 500 Oe, we estimate for a 50 A x 1/xm domain 
cross section a flux of approximately one flux quantum. 
In order to average the influence of domains, we have 
measured different junctions with the same thickness of 
the ferromagnet layer. 

Fig. 3 shows the variation of Tc in the Ni-Nb-Ni(x)- 
Nb-Ni multilayer structure versus the thickness of the 
Ni(x) layer. The data contains a pronounced minimum 
around Xm\n—^'^ A, which is in excellent agreement with 
the expected ^ tt transition at dmin ~15-17 A implied 
by the critical current measurements above. In order to 
ensure that the observed minimum arises from a — > tt 
transition, the variation of Tc in a Ni-Nb-Ni(x) structure 
was measured for the same range of Ni(x) thicknesses, see 
Fig. 3. For symmetry reasons, Tc should vary twice as 
fast for the trilayer compared to the 0-phase of the pen- 
talayer. The absence of a pronounced local minimum in 
the trilayer with a single Nb layer therefore undoubtedly 
indicates that the minimum observed for the pentalayer 
containing two Nb layers must arise from a ^ tt tran- 
sition. The kink- like change of Tc at 17 A also suggest 
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FIG. 3: Critical temperature Tc as a function of Ni-layer 
thickness xni for the trilayer and pentalayer structures. The 
experimental data are shown as symbols, while the curves are 
fits to the data with the theory discussed in the text. For the 
pentalayer, the highest Tc is obtained for a zero junction (full 
line) for dMi <16 A, and for a 7r-junction (dashed line) in the 
range 16 A< dwi <40 A. 

that there is a crossing point of two phases, namely the 
0— and TT-phases. 

In order to compare our experimental results for Tc 
with theory we solve the gap equation and compute Tc 
of our systems with the quasiclassical Green's function 
technique in the diffusive approximation 21] . Near Tc, 
the order parameter A <C Tc and the Usadel equation 
can be linearized. We have generalized the results for 
symmetric trilayers by Fominov tt al. in Ref. [2^ to 
asymmetric F1-S-F2 trilayers and symmetric F2-S-F1-S- 
F2 pentalayers. Instead of discretizing the spatial coor- 
dinate we Fourier-series expand the order parameter and 
flnd Tc by studying the resulting eigenvalues of the gap 
equation. With this technique |23l |. the accuracy as well 
as the speed of the numerics are improved immensely 
compared with previously used methods |2^ 01 • 

Several material parameters serve as input to the 
model: the exchange fleld Ecx of Ni, the critical tem- 
perature Tco of Nb in the absence of the Ni layers, and 
the diffusion constants of Ni (Dp) and Nb (T'g). The 
boundary conditions [2^ at the Ni-Nb interface are ex- 
pressed in terms of a normalized boundary resistance 
and the conductivity mismatch 7 between the Ni and Nb 
materials. We have considered these two quantities as 
free parameters. 

Our flts of the experimental data for Tc as function 
of Ni layer thickness with the above theory are shown 
as curves in Fig. 3. We use as input parameters the ex- 
change fleld Ecx = 200 meV and the diffusion constants 
Dp = 2.8 cm^/s and = 3.9 cm^/s (with D = ^Vfi) 
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obtained from the fit of Ic in Fig. 2. The fit parameters 
are the bulk Nb transition temperature Tco, the interface 
resistance 7;,, and the materials' conductivity mismatch 
7. The fits indicate that T^o of the batches of trilayers 
and pentalayers differ, while other sample characteris- 
tics remained essentially the same. Although all samples 
within the trilayer set and within the pentalayer set were 
prepared in situ, both sets were evaporated separately. 
We assign the difference in Tco to this fact. 

We have calculated Tc as a function of x^; for the pen- 
talayer for zero phase difference and for tt phase difference 
between the two superconductors, using the same param- 
eters. The corresponding curves are shown in Fig. 3 as 
full and dashed lines. The ^ tt transition takes place 
where the two Tc curves cross. We note, that the fit pa- 
rameter Tco is determined by the small thickness data 
points, and the remaining fit parameters 7 and 76 are 
determined by the fitting of the zero-phase curve. Hav- 
ing no additional fit parameter, the — > tt crossing at 
XNi = 16 A is in remarkable agreement with the experi- 
mental data for Tc, and with the prediction of the Ic data 
fit. We also fit the experimental data for the trilayer with 
very similar interface parameters. 

We would like to mention that we expect corrections 
to the Usadel theory when the exchange field is large. 
From the fit of the Ic data, where such corrections were 
taken into account, we see that £ ~ while Usadel 
theory works well for £ <C Lm- Nevertheless, we obtain 
a remarkably good fit for Tc as function of xn;. This is 
probably due to the fact that the Ni-layers for the Tc 
measurements are quite thin, in which case surface dis- 
order is relevant and justifies the use of Usadel theory. 
We note that the surfaces are characterized in this case 
by strong disorder with a large number of point contacts 
(high-transmission channels) . This is consistent with our 
observation of regions with short-cuts in the samples used 
for the Tc measurements, which prevented us from ex- 
tending our Ic measurements to d^; < 30 A. 

In summary we have demonstrated that both the crit- 
ical Josephson current and the critical temperature of 
Nb-Ni multilayers vary with the Ni thickness with ap- 
proximately the same period, 16±1 A. We deduce from 
the period a magnetic length Lm=^0 A, corresponding 
to an exchange energy of Ecx='200 meV. By measuring 
Tc in Ni-Nb pentalayers, we have observed a ^ tt tran- 
sition at a thickness for the central Ni layer consistent 
with the theoretical prediction using Usadel theory. For 
higher thicknesses, we see further tt — )■ and ^ tt tran- 
sitions in the critical Josephson current, consistent with 
the period in our Tc measurements and with the predic- 
tions of theory. Our results demonstrate the feasibility 
of using strong ferromagnetic materials in the design of 
Josephson devices for future applications. 
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